 1:


 a := 0; 
 evalf(a); 
 b := 0; 
 c := 1; 
 E := a*z^2+b*z+c; 
 niy := .4; 
 G := E/(2*(1+niy)); 
 Alfa11 := 1/E; 
 Alfa12 := -(1*niy)/E; 
 Alfa13 := -(1*niy)/E; 
 Alfa22 := 1/E; 
 Alfa23 := -(1*niy)/E; 
 Alfa33 := 1/E; 
 Alfa44 := 1/G; 
 Alfa55 := 1/G; 
 Alfa66 := 1/G; 
 
 Liambda := 2*niy*G/(1-2*niy); 
 Liambda11 := 2*G+Liambda; 
 Liambda12 := Liambda; 
 Liambda13 := Liambda; 
 Liambda22 := 2*G+Liambda; 
 Liambda23 := Liambda; 
 Liambda33 := 2*G+Liambda; 
 Liambda44 := G; 
 Liambda55 := G; 
 Liambda66 := G; 
 
 A := Matrix(3, 3, [[Alfa11, Alfa12, Alfa13], [Alfa12, Alfa22, Alfa23], [Alfa13, Alfa23, Alfa33]]); 
 B := evalm(1/A); 
 
 Xalfa := -1/G; 
 Yalfa := -1/G; 
 Zalfa := -1/(2*G+Liambda); 
 
 a1 := -(2*G+Liambda)/G; 
 b1 := -1; 
 c1 := -(G+Liambda)/G; 
 d1 := -1-Liambda/G; 
 
 a2 := -1; 
 b2 := -(2*G+Liambda)/G; 
 c2 := -(G+Liambda)/G; 
 d2 := -1-Liambda/G; 
 
 a3 := -(G+Liambda)/(2*G+Liambda); 
 b3 := -(G+Liambda)/(2*G+Liambda); 
 c3 := -G/(2*G+Liambda); 
 d3 := -G/(2*G+Liambda);

 2:


 beta13 := B[1, 3];
 beta23 := B[2, 3]; beta33 := B[3, 3];
 z := 0; evalf(beta13); evalf(beta23); evalf(beta33); evalf(Liambda44); evalf(Liambda55);
 z := 0.1; evalf(beta13); evalf(beta23); evalf(beta33); evalf(Liambda44); evalf(Liambda55);



